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Abstract 

A dynamical model of the natural conflict triad is investigated. The 
O I conflict interacting substances of the triad are: some biological popula- 

(N ■ tion, a living resource, and a negative factor (e.g., infection diseases). We 

^ ■ suppose that each substance is multi-component. The main coexistence 

phases for substances are established: the equilibrium point (stable state), 
the local cyclic orbits (attractors), the global periodic oscillating trajecto- 
ries, and the evolution close to chaotic. The bifurcation points and obvious 
thresholds between phases are exhibited in the computer simulations. 
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g : 1 Introduction 



We use the term conflict triad for notation of a physical system consisting of 
three conflict substances (denoted by P, R, Q) which exist in the same space fl. 
Each substance at initial moment of time is presented by a real value P, R, or 
5^ ! Q in accordance with its amount characteristic. The whole system is complex 
since every substance contains a family of components distributed along regions 
Qi which compose the existence space: Q = IJILi 2 < n < oo. So one can 
think that P, R, Q are vectors with non- negative coordinates: 

P = (Pi,...,P„), R= (i?i,...,P„), Q = (Qi,...,Q„), Pi,R^,Q^>0, z = l,...,n, 

P = Pl + --- + Pn, R = Rl + --- + Rn, Q = Ql + --- + Qn. 

We study the evolution of the system at discrete time: 

{P^, R^, Q^} ^ {P^+\ R^+\ Q^+i}, iV = 0, 1..., 

where P*^ = P, HP = R, Q° = Q and the map is defined by the law of conflict 
interaction. In general this law is unknown. Our definition of the transforma- 
tion * (see below formulae fl3.17p . (13.181) ) starts with the famous Lotka-Volterra 
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predator-prey approach (see e.g. (TUl [HI [221 [21]) and is based on heuristic un- 
derstanding of the physical essence of the substances P, R, Q. 

The schematic picture of mutual dependencies between P, R, Q may be pre- 
sented by the following diagram: 




where each signed arrow means a direction of positive or negative influence one 
of substances at another. The positive dependence involves a certain growth and 
negative - decreasing of quantitative characteristic of the corresponding specie. 

At this diagram we have three types of interaction connections: two variants 
of the "plus-minus" interactions and a single version of the "minus-minus" 
interdependence . 

The "plus-minus" type of interaction is in fact the modified stochastic version 
of the predator - prey model (c.f. with |25]). Here we use two different variants of 
this type. First, the conflict interdependence between behavior of some biological 
population playing a role of predator and the dynamical changes of the vital 
resource as a prey. Second, the nonlinear relations between the local density of 
biological species (as some kind of a prey) and the expansion processes of viral 
infections, like a certain predator. 

The "minus-minus" interaction is essentially different. This type describes 
the conflict competition of mutually alternative substances like virus infection 
diseases and medicine. In the paper we consider this couple slightly wider, we 
treat the vital resources as the opposite substance to the viral infection and admit 
its quantitative decreasing under destroying actions of the latter. 

In more details, one can understand a sequence Pj^, = 1, 2, ... as the quan- 
titative evolution of some biological specie in a fixed region Qi under influence 
of two opposite factors. Namely, the population dynamics of Pf^ is determined 
by the positive dependence produced by the vital resource environment from the 
side R that provides some growth of Pf^ , and by the negative influence gener- 
ated by Q. The possible viral infection caused by the substance Q leads to the 
quantitative losses for Pf. Thus, at the moment + 1 the value P^~^^ has to 
be proportional to amount of the biological spices in region f2j at the previous 
moment, it should grow at increasing of Rf, but to fall at increasing of Qf. 

In turn, under assumption the vital resource is not exhaustible globally, the 
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regional changes of its amount Rf are dependent naturally from the local (=re- 
gional) intensity of their using as a source for existence of biological species, and 
in addition from the negative influences of Q which has the alternative nature. 
Finally, the evolution of coordinates Qf describes the dynamical picture for be- 
havior of threat concentrations for existence of biological species inside every 
separate region So the local expansion of some infection Q^'^^ grows at 
increasing of Pf^ and decreasing of Rf^ . 

The explicit formulae of all mutual interdependence will be defined in Section 

3. 

One can associate the triple of substances P, R, Q with the philosophical 
triad: mankind, good, and evil or, in accordance with mythology, with flora 
and fauna of Earth surrounded by water and fire as a positive and negative 
resources. More specifically, may mark the quantitative global population 
of some biological specie (humanity), R^ - the variation of vital resources in 
the current environment, and corresponds to the dynamics of threatening 
concentrations associated with various negative phenomena that may influence 
both at the biological population P^ and the resource environment R^ too. In 
particular, coordinates i^/^, Rf, Qf describe the local evolution of components of 
triple indicated substances in the fixed region VLi. 

The general scheme of the dynamical system studied here is complex enough 
and it is hardly to expect for establishment of abstract mathematical results 
(theorems) concerning the global behavior. Nevertheless, we got some exact 
results which arc true, in particular, for separate two-sided links inside of the 
conflict triad, i.e., for couple interaction between flxed substances (see Section 2). 
We study the whole complex system by using the series of computer models in 
Section 3. Their analysis exhibits remarkable features which would be useful in 
applications. 

Two intrinsic questions arise. Whether the construction of conflict triad is 
well-defined? If so, whether the dynamical picture of the substances P, R, Q 
coincides with our intuitive pattern of the conflict triad coming from natural 
practical experience? 

The flrst important result of this work is the setting of mathematical consis- 
tency of our constructions. All the interrelations used in the dynamical system of 
the conflict triad do not lead to collapse. More exactly, it is shown that under a 
certain choice of additional parameters of the model (this choice is the non-trivial 
problem itself) the whole complex system does not decay (destroy themselves). 
This means that the values of all coordinates are changed in the physically rea- 
sonable scopes. 

The second important result of our work is that various computer models 
that realize the conflict triad expose the typical features for the complex systems 
behavior. Namely, analysis of the values P/^, Rf , Qf in concrete models demon- 
strate the distinctive properties that are observed in natural conditions where 
more than two substances take part in confrontations. Among such properties, in 
particular, we observe the presence of distinct phases for existence of the system: 
the phase of dynamic equilibrium (the equilibrium state), the phase of periodic 
oscillations which contains the wide spectrum of trajectories that have tendency 
to approximate the cyclic attr actors, the presence of areas for bifurcation points 
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and thresholds between different phases, finally we meet the class of evolution 
close to chaotic (which is called the quasi-chaotic). 

We note that mathematical basis of our constructions is a concept of the 
conflict dynamical system that has been developed in [121 1131 El [13 121 El US] , 
(c.f. with Among other deep theoretical researches that influenced on our 

work we refer P S] , [IB] - [25] . 

The most close to our conflict triad is the well-known SIR-model and its 
different variants that describe the dynamic of epidemic infections (for details 
see, for example, [221 123]). In the SIR-model the idea of conflict triad is actually 
presented but in an implicit form. So, all the population of some biological 
specie is divided into three groups: S is the amount of persons favorable to the 
infection (susceptible), / - infected and able to carry the infection, and R - which 
are under cover or recovered already. The complicated SIR-model contains the 
additional group E of those, who has an infection in the hidden (latent) form. 
The evolution of the above mentioned groups in time is determined by relatively 
simple but nonlinear equations: 

S = —rSI, I = rSI — al, R = al, 

where coefficients r, a > 0. It is assumed that S(t) + I(t) + R(t) = const and since 
that 5" -|- / + i? = 0. The main problem is to research the dynamic of distribution 
for the infection I{t), in particular, when I{t) increases or decreases depending 
on coefficients r, a and starting values S'(O), /(O), i?(0). There exist a number 
of publications (see references in [221 l23l E]) where the SIR-model is improved 
or modified in accordance with concrete specifics of biological species, view of 
epidemic infections, and conditions of its expansion in a certain environment. 

In comparison with the SIR-model, our construction of the conflict triad is 
substantially more perfect in two principal aspects. 

First, we use a partition of whole existence space into family of finitely many 
regions. Of course, it corresponds to the widely observed environment phe- 
nomenon of nature: the existence space always is separated into bounded do- 
mains. That is why the computer models of the system are able, in particular, 
to describe the expansion of disease infections of biological population separately 
in every region of existence and to define, for example, the most safe regions. 

Another new important step in our work is that we put as a basic the prob- 
abilistic (statistical) law of interaction between conflict substances. We consider 
that statistical interdependence reflect more deep connections between opponents 
any nature. Thus we may describe not only quantitative changes of populations 
but also its redistribution along regions of existence caused by the conflict inter- 
action. 

Shortly, in contrast to the above mentioned SIR-model our method of the 
construction of the conflict triad uses a partition of the existence space into 
regions and application of statistical formulas for description of the dynamical 
picture. 

Finally, it is especially important, we assume that each substance (opponent) 
is a priori non-annihilating. This principle is putting directly in formulas f l3.17p . 
03.181) that govern the conflict interactions. So, any of substances P, R, Q can 
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not annihilate another one, but has to reach the compromise equihbrium state or 
migrate along regions by some law. 

2 Dynamics of the bilateral couple interactions 

According to equations (I3.17p . (13.181) (see below) the interactions between sub- 
stances P, R, Q in the conflict triad are nonlinear and rather complex. But the 
separate bilateral couple interdependency admits rigorous enough mathematical 
analysis. In this section we state two results for situations marked as "plus - 
minus" and "minus - minus" models. 

In the first bilateral "plus - minus" model we analyze the dynamics of changes 
for substances P and R. Plus means that P is positively influenced by R, and 
minus - that R is negatively influenced by P. That is a typical predator - prey 
situation (see for example [21]). However we consider different, more specific 
functional dependence between opponents. 

Let P^, = 0, 1, denote the amount evolution of P (some fixed biological 
specie) at discrete time. As soon as the space fl is divided onto family of separate 
regions Q = IJiLi^*' 2 < n < oo, the complete amount at every time 
moment is certainly distributed along these regions: = P/^ + ■ ■ ■ + ■ For 
simplification of the problem, we assume that the complete amount of substance 
at the whole territory Q is permanent P^ = P = const. It means that the 
average rate of growth and the decay rate of the substance P population (the 
rates of birth and deaths) are independent by the time, i.e. are in the dynamical 
balance. In fact this occurs often enough, at least locally at certain periods of 
development of a complex physical system. Thus, in the such simplified model 
we analyze only the redistribution dynamic for Pf^ along the regions fli. In 
the model of conflict triad the redistribution of values Pf^ is stipulated by two 
factors: the positive influence of the substance R (the population of Pf^ grows 
using the vital resource environment), and the negative influence from the side 
of Q (for example, threats, infection diseases cause some decay). In a situation, 
where the global influence of the last factor is insignificant and negligible, the 
bilateral coupled interactions between regional values of P^ , R'^ are essentially 
simplified in the mathematical sense. In turn, the dependence of R from P is 
purely negative (the amount of the substance R is "burned" as a vital resource for 
P). However we assume the global amount of P in f2 is stable R^ = R = const. 
This stability condition means that R is continuously supplied due to the external 
source (such as the Sun) that now is not examined. Nevertheless at each time 
moment the regional distribution of the vital resource R^ = R^ + ■ ■ ■ + R^ is 
changed according to the appropriate law. 

Thus the simplified problem is to study the redistribution dynamic for values 
Pj^, Pf' along regions f2j produced by a certain "plus - minus" interaction 
between substances P and R under assumption of the global amount stability for 
P, R in Q and their independence with the third substance. 

The proper conflict dynamical system has the view: 

{pr, pf } ^ {pr\ R^^'}, = p^. = R^- (2.1) 
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Surely it is considerably simpler as compared to the behavior of all complex 
system. After the explicit definition of the conflict transformation * (see below 
(I2.3P ) we are able to fulfill the detailed enough analysis of system (12.11) . The 
main results are stated in Theorem 1 which gives complete enough description of 
behavior of this system. 

Let us introduce the simplest variant of concrete formulas of the "plus - 
minus" interaction between substances P and R . We write down these formulas 
in terms of coordinates of stochastic vectors with a unite /i-norm: 

pN tdN 

p^ = (pf...pf,...,p:), r^ = (rf...,rf,...,r:), pf := i^, rf := ^, (2.2) 

where, we remind, P, R are the amount characteristics of substances P, R in 
the whole space fi. Namely, at + 1 step the coordinates p^~^^ and r^~^^ are 
determined iteratively by the rule 

pf+i = ^^il±^, rf-^^ = ^^^^^, N = 0,1... (2.3) 

where the normative denominators 

<. = l + (P^r^)' 4 = i-(p^r^) 

((■,■) denotes the inner product in R"') ensure that vectors p^"''^,r^+^ remain 
stochastic. 

We note that different signs in the numerators of formulas (12. Sp just determine 
the essence of the "plus - minus" model. Plus means that P^^^ increases 
depending on a value Rf , and minus supplies the decay of Rf^^ at increasing Pl^ . 
Of course, one have to make these interpretations after re-normalizing inverse to 
(12. 2p . In fact, the model is well-defined by virtue of the normative denominators. 



Theorem 2.1 ("Plus-minus" model) The conflict dynamical system ^2l\) given 
by the formulae ^2.3) has three typical phases of behavior. 

The first phase (equilibrium point) is determined by the uniform distribution: 

yN Pl^ = P/n, Rf = R/n- p = (1/n..., 1/n), r = (1/n..., l/n). (2.4) 

However this equilibrium state is unstable. 

The second phase (existence of the limiting fixed points) is stipulated by the 
condition: at least for one i, 

Pi = r^^O. (2.5) 
In this case the system trajectory converges to the limiting stable state, 

p°° = lim p^, r°° = lim r^, p°° ± r~, (2.6) 

which is invariant with respect to the conflict interaction. 

The third typical phase (the quasi-chaotic behavior) occurs under the starting 
conditions 

p^r, Vz,p, T^Or, 7^0. (2.7) 

At this phase each of coordinates pf ,rf i = l...,n oscillates between zero and 
one, in general case without any regular law. 
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Proof. In the case of uniform starting distributions of P, R along regions, 
Pi = P/n, Ri = R/n,yi, all coordinates of the stochastic vectors p, r are equal: 
Pi = Ti = 1/n. Then it is easy to find that 

(p,r) = l/n = Zp^r = Zr,p- 

Therefore due to (12. 3 p for all we have pf = rf = 1/n. It proves that uniform 
starting distributions define the equilibrium point for dynamical system ( 12. ip . 
This state is unstable. An arbitrary small deviation e > of any coordinate Pi 
or Tj from 1/n leads in time to greater deviations (see Lemma 1 below). 
Let us prove (12. 6 p under condition (12. 5p . Here we introduce the value 

i 

and call it the conflict index for dynamical system at moment N. 

Without loss of generality we assume that (12. 5 p is fulfilled for only single co- 
ordinate and all other ones, ri,pkjki, r^^j are nonzero. Then, it follows from (12. 3p . 
we have to show that there exist the limiting vectors p°°, r°° with coordinates: 

pr = 0,C = l, p->0,r~ = 0, ky^i. 

Indeed, if = and ^ 0, then due to < 6^ = < 1 the sequence rf = 
r^~^/(l — 9^^^), N = 1,2,... monotonically increases. Thus, since < 1, 
there exists the limit r°° = limAr^oo "'"f^- That is, because = ■ Y[h=i V(l ~ 
6'') < 1, the convergence of rf' implies with necessity that the conflict index 6'^ 
monotonically decreases. In fact — ?• 0. To see rf° = 1 we prove that all 
r^-^- — )■ 0, N ^ oo. Assume the opposite, i.e., that there exists at least single 
coordinate , k ^ i which does not converge to zero. Then, due to the conflict 
index 6^ ^ 0, the coordinate p^ have to come to zero. By (12. 3 p we have 

Pr^=pf ■(l + rf)/(l + ^^^)^0, 

that is possible if only < 9^ , thus it is the contradiction. So — )• 0. And 
by the same reason all other coordinates converge to zero too. Therefore 

= 1. By similar way we obtain the existence of the limiting coordinates 
p^, k ^ i, which may take any non-zero values given unit in a sum. It is evident 
due to 6°° = that p°° _L r°° and therefore in the limiting state the system 
reaches the stable equilibrium. 

Let us consider the case of the third phase. We assume that vectors p, r are 
different and all starting coordinates are non-zero. Then < ^ < 1, the vectors 
p, r are not orthogonal. In particular, we exclude that Pi = = 1 for some i. Let 
us show that in such a case all coordinates oscillate inside open interval between 
zero and unit (without any evident pattern). 

Our argumentation is based entirely on formulae (12. 3p . Take any couple of 
coordinates ^ pi ^ Vi ^ 0. Assume that in the starting moment the inequalities 

ri<p^<e (2.8) 
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hold. We shall show that all other possible inequalities between r^,pf,6^ will 
successively appear at some moments of time. In the next we use the evident fact 
about the qualitative behavior of values 9^, pf ■ Namely, by the definition 
of the conflict index, = '^^pf ■ rf , its small changes are slower than changes 
of any fixed coordinate pf,rf. Indeed, a finite sum of the differential products 
ApjAr, has the second power of smallness with respect to any separate differential 
Api or Arj. In what follows we omit the time subscript A^. 

Directly from (12. 8p by (12. 3 p it follows that pi decreases and rj increases. This 
leads with necessity to the transformation (12. Sp into the inequahty 



Pi <r,<9. (2.9) 

Using again formulae (12. 3 p we check that the coordinate Pi will be else decreasing, 
and Ti - increasing. This is continued till the moment when instead (12. 9p appears 
the inequality 

Pi <9 < Ti. (2.10) 

In turn, again by formulae (12. 3p the coordinate pi begins to grow due to (I2.10p . 
although Tj is still increasing. This tendency continue till the moment when the 
inequalities (l2.1Up change at 

9<p.,<n. (2.11) 

We note, inequalities (12. lip do not imply the convergent to zero of the conflict 
index, i.e., the vectors p, r could not become orthogonal. 

By (I2.10p the coordinate Pi grows. It continues to grow after coming to (12. lip , 
but Tj begins to decay. This leads to the new inequalities 

9<ri< Pi. (2.12) 

In turn, it follows that pi still grows, and continues to decay. On this way the 
inequalities 

r^<9<p^ (2.13) 

appear. The latter produce the decaying of Pi and also decreases as at previous 
period. But these changes continue only till the moment when (I2.13P is replaced 
by the starting inequahties (12. 8p . 

We note that the sequence of transformations from (12. Sp to (12.130 are ordered 
and sometimes equalities may appear but by (12. 3p they pass into inequalities 
immediately at the next moment. Thus the full cycle of all possible inequalities 
between 9,ri,pi was realized. 

By the way we observe that no one coordinate ri,pi, as well as the conflict 
index 9 could not go to zero closely. In particular, for example, if (12. 8 p takes place 
then Ti begins to grow. In general, as soon as some coordinate becomes smallest 
it begins to increase with necessity that follows from (12.31) . 

It is not hard to see that for different couples rj,pj each inequality from (12. 8p 
- (I2.13P is fulfilled at various moments of time and the corresponding transforma- 
tions are not synchronous. That is why the existence of regular cyclic oscillations 
for the all complex system in the considered situation is questionable. Therefore 
in this phase orbits of the system are similar to some kind of the quasi-chaotic 
behavior. 
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If one or several coordinates = ■ ■ ■ = = 0, 1 < m < ra, but all Pi ^ 
0, i = 1, ...,n, then it is easy to see that p^, ■■■,Pi^ 0, and residual coordinates 
have the quasi-chaotic behavior. 

The proof will be completed if we prove that the equilibrium state (12. 4p is not 
a one-point attractor (see Lemma [2.21 bellow). 
□ 

Lemma 2.2 An arbitrary small deviation from the equilibrium state ( [^.^[ j of the 
conflict dynamical system h2.1\) : 

Pi = l/n ^ p^e = ^/n + Ei, rj = 1/n r 5 = 1/n + ^£^ = 0, y^(5» = 

i i 

automatically leads to the greater deviation at least for some coordinates: 

Ei-f e\> Ei 5i 5i > 5i. 

Proof is not trivial and here we present only its sketch. We need to analyze the 
dependence of values pf = rf from fixed deviations as — )• oo. The linearization 
of (12.31) shows that the terms of the first order by Ei,5i expose the following 
changes of deviations 

Ei-^ E[ = Ei + 1/ (n + l)5i, 5i -> 5[ = 5i-ll{n- l)Ei. 

These formulae are defined by the strictly positive definite matrix t with elements 
^11 = ^22 = 1, ti2 = ^/{n — 1), ^21 = —^/{n — 1) which does not depend on the 
starting deviation. By this reason the iteration of formulas (12.31) produce new 
deviations increasing with N ^ oo, that one may check directly. 
□ 

Computer simulations also confirms that the equilibrium point is not stable. 

In a general situation the above considerations exhibit the existence of infinite 
oscillations for non-zero coordinates. Whether these oscillations may be cyclic? 
That is, whether the cycles have a finite number of steps? Apparently it is 
possible under a certain starting connection between values of all coordinates. 
However, the existence of exact finite cycles in the "plus - minus" model is the 
open question until now. 

The similar characteristic behavior has the conflict dynamical system with 
the bilateral interaction of "minus - plus" type - the model which describes 
the evolution of biological species under influence of some infection (the viral 
environment). 

In such a case we have to define a vector of an initial statistical distribution 
for a virus infection along regions flf. 

q = (gi...,gi, ...,g„), qi:=—, 
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where Q = Qi + ■ ■ ■ + Qn- Then the evolution changes of pf, are governed by 
the formulae similar to (12. 3p : 

Pi — M ) % - AT ) l^-J-^J 

where normalizing denominators z^^ = 1 — (p^, q^), z^^ = 1 + (q^, p^) ensure 
that vectors p^, are stochastic. The opposite signs in the numerators of fl2.14p 
have now the following interpretation: minus means the decreasing of a biological 
population caused by infection, plus provides the growth of virus concentrations 
in regions with large local amount of biological species. In fact, according to 
(I2.14p . the full amount of biological species and average virus concentrations are 
stable because normative denominators z^^, z^^ guarantee the non-annihilation 
of a population and the natural dissipation for bacteria. Of course, the real 
quantitative changes of these substances in various regions are determined by 
fl3.17p under the complex triple conflict interaction. 

Let us now consider an abstract variant of the "minus - minus" model. It 
may be interpreted as a situation of conflict fighting between couple of purely 
alternative opponents of type "infection - medicine". Shortly this alternative 
confrontation may be written as "either - or" that describes a tendency to 
exclusion one to other from every region. In the terms of stochastic vectors r, q 
the alternative interaction we represent by the following formulae: 



N+i _ rf) j^j^^ _ rf{l qf) 



where a value of the normalizing denominator z'^ = 1 — 6^ strongly dependents 
of the conflict index 6^ = (q^, r^). 



Theorem 2.3 ("Minus - minus" model) Given a couple of stochastic vectors 
q, r G M",n > 1 assume the conflict index 9 = (q, r) satisfies the inequalities: 

0<9 <1. 

Then each trajectory of the conflict dynamical system 

{q^, r^} ^ {q^+\ r^+^}, q° = q, r° = r, iV = 0, 1... 
generated by Ii2.15\) goes with necessity to the equilibrium state: 

q°° = lim q^, r°° = lim 

which is a fixed point: q°° = q°°*r°°, t°° = r°°*q°°. Moreover, the limiting 
vectors are orthogonal q°° ± r°°, if q 7^ r, and identical, q°° = r°° , if the 
initial vectors are equal q = r. In latter case the coordinates of limiting vectors 
are uniformly distributed: q°° = r°° = 1/m, where m < n denotes the amount of 
non-zero initial coordinates. 
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For the proof see [HI [13] and [H US El El [7] . 

We remark that in [2] the formulae fl2.15p was used for the construction of 
complex system which generalize the well-known predator-prey model and is in 
fact some variant of the vector analog of Lotka-Volterra equations. Here we 
mention the paper [5] where the idea of clusters (districts, regions) was also used 
in three-dimensional discrete-time Lotka-Volterra models. 

Recall, that here we construct the discrete time models. However, just below 
we exhibit a pair formulae with continuous time: 

./ X p(0 + r) , , r(9 — p) 

where p = p{x,t),r = r{x,t), x ^ Q denote the distribution densities of the 
corresponding substances and nip^r = P ■ R, 0{t) = J^p{x,t)r{x,t)dx. We plan 
to study these equations in consequent publications. 



3 The conflict triad 

We write the conflict dynamical system that describes the evolution of simulta- 
neously interacting triple substances P, R, Q as follows: 

{P^,R^Q^}^{P^+\R^+\g^+i}, N = 0,1... (3.16) 

where the conflict map (composition) * is defined below by formulae f l3.17p . 
(13.181) . As above we assume that substances P,R, Q have a common space of 
existence fl, which is decomposed in a natural way into the finite set of separate 
regions, Q = IJiLi^*' > 2. The conflict triad is a complex system. It means 
that each substance has an inner structure: 

P = (Pi...,P„), R= (i?i...,i?„), Q = (Qi...,Qn) 

where the elements Pi,Ri,Qi i = 1., ...,n determine the proper quantitative de- 
scription of the corresponding substances. 

The investigated substances of conflict triad have different physical nature. 
That is why the concrete formulae of interactions of every substance with a com- 
plementary pair, namely P with the pair {R, Q}, R with {P, Q}, and Q with 
{P, R} are essentially different one from another. We represent the complete 
mechanism of interconnection that is contained in the conflict composition * into 
two parts: formulae (13.171) which gives the algorithm of quantitative changes of 
absolute values Pi,Ri, Qi in regions fli and f l3.18p which describes the statistical 
law of redistribution of occupation probabilities of regions Qi by substances 
P,R, Q. 

The evolution of quantitative regional changes for Pi,Ri, Qi is assigned by the 
equations: 



R 
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= ^'^"^'fl = P.R^ = R.Q',=Q.N = 0, 1..., 

where parameters di^d^id^ > characterize the rate of changes intrinsic to the 
real model. The normalizing denominators Zp , Zq , Z^ ensure the stable global 
amount of the proper substance in the whole space fl. Of course, the global 
amount of every substance may be additionally changed by virtue of external 
circumstances, but we do not consider such influence. Here, for the sake of sim- 
plicity we assume that the total quantitative characteristics for each of substances 
are unchanged and therefore we may write: 

n n n 

P = 5^P. = P^ Q = $^Qi = Q^, /? = ^i?. = /?^, N = o,i... 

i i i 

We interpret formulae f l3.17p as follows. The quantitative growth of the bi- 
ological species in Qi region on + 1 step of the conflict fight is proportional 
to amount Pf^ at the previous moment of time and to the difference values with 
some coefficient of the vital resource Rf and the factor of elimination threats Qf. 
At some moment of time N it may happen that above difference has negative 
value. Then Pf^ will decrease quickly enough. However, such period of devel- 
opment has to be short. Otherwise, the system will be destroyed and loses its 
physical sense, for example, if some of coordinates Pf^ , Qf , Rf becomes nega- 
tive. Similarly, we interpret the dependence of Qf from the same difference, but 
with another coefficient. In turn the quantitative changes of the vital resource 
Rf are very sensitive to the relative density of threat for existence of biological 
population Qf / P^^ . In real models the coefficient is small. We note that for- 
mally, according to numerators in formulae fl3.17p all coordinates are increasing. 
Nevertheless due to the normalizing denominators the dissipation process courses. 
This automatically provides the decreasing of all values P^^, Rf , Qf at each step 
of the conflict fight. Of course, the concrete character of interdependencies, their 
physical interpretation, and the role of parameters is determined by the model of 
research. 

The second part of our mechanism of the conflict interaction * has purely 
probabilistic inter-regional character. To write it in the mathematical terms it is 
necessary to transfer the vectors P, R, Q into stochastic ones: 

p = r = (ri...,r„), q = (gi-.., g^), 

where the coordinates 

Pi := Pi/P, Ti = Ri/R qi := Qi/Q i = l...,n 

have a sense of probabilities to flnd the corresponding substance P, R or Q in i-th 
region. In other words they are the occupation probabilities of fli by P, R, Q. 
In the course of interactions the redistribution of these probabilities takes place. 
The law of these changes is determined by the following formulae. They are some 
statistical variants of Lotka-Volterra discrete time equations flUi [TT] (c.f. with 

1): 



Pfa + airf-qf)) 
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where parameters a, 5, c > characterize the intensity of the conflict redistribu- 
tion. 

So the first of these formulae shows that the statistical redistribution for the 
population substance P is maximal in the region with the highest probability 
to find the vital resource and the lowest threat to existence. In turn the sec- 
ond formula implies the decreasing of a probability to find the vital resource 
R in a region where the biological population is large (because the later uti- 
lizes the vital resource) and there is a high statistical infection concentration 
Q. Finally, according to the third formula in f l3.18p the probability of the in- 
fection threat (for the population existence in i-th region) increases together 
with growing of population and decreases under action of the vital resource as 
an alternative substance. The denominators in (13.181) provide that all vectors 
= (pf,...,p^), = «,...,r^), = (gf,...,g^), iV = 1, 2, ... have 

unite norms. 

To complete the definition of the conflict map it is necessary to fulfill the 
re-normalizing of the vectors p^"*"^, r^+^, q^"*"^ after using the formulae (13.171) 
and (13.181) . It means transferring to the final quantitative values of substances 
P, Q and R at each {N + l)-th step: 

pr' = p^' ■ P. Q^' = ■ Q, = ■ r- (3.19) 

Thus, the conflict composition as a map * in (I3.16P that generates the dy- 
namical system of conflict triad is entirely determined by formulas (13 . 1 71) . (I3.18p . 

In the present work we made in fact only the first attempt to construct and 
analyze the simplest computer models of the conflict triad. However even this 
activity finds out the series of interesting observations usually inherent to com- 
plex systems. In particular, we establish the existence of fixed points (which are 
attractors), the existence of the stable limiting equilibrium states, the appear- 
ance of cyclic orbits, which are attractors too, the critical bifurcation points, the 
oscillating trajectories without an evident law of behavior, most similar to quasi- 
chaotic. So, we hope, the subsequent research will lead to series of more deep 
results and useful applications. 

3.1 Computer models. The equilibrium state 

Due to formulae (I3.17p . (13.180 . (I3.19P the conflict triad is the highly complex 
system. That is why the question of existence of its state in which the appropriate 
substances are situated in the equilibrium balance does not have an obvious 
answer. In physical reality we usually observe that complex systems possess such 
states as a rule. We mean the existence of the dynamical equilibrium in a complex 
picture of the interaction process between opponents with contradict tendencies 
when the various alternative substances coexist. We remark that the presence 



13 




Figure 1: The equilibrium fixed point in the phase space {Pf, Qi) 

of such state does not follow straightly from Theorems 3.1 and 3.2 ensuring the 
conditions for existence of the equilibrium states. However one can conjecture 
about the possible mutual compensation of oscillations which inherent to the 
models with bilateral "plus- minus" interaction. It has to lead to stabilization 
and to the equilibrium state as soon as the limiting compromise distributions 
exist. Some analogs of such state are found in Theorem 3.2 for "minus-minus" 
models. 



Theorem 3.1 The conflict triad dynamical system l[3.16]) defined by Iji3.n\ ) l[3.18\) 



^3.iy\) . under the condition that all starting coordinates are non-zero, possesses 
the equilibrium state. This state is determined by the fixed point with coordinates 
P^'^,R^^,Q'^'^,i = l...,n which are equal to the arithmetic mean values of starting 
amounts of substances P, R, Q in regions Qi. 

Proof. By virtue the assumption that global amounts P, R, Q are constant, the 
coordinates of any vector P, R, Q cannot increase or decrease simultaneously. 
So, if we assume that at least one of coordinate is changed, for example, becomes 
bigger, then there exists another one which will decrease with necessity. However 
it is impossible since by fl3.17p . f l3.18p . fl3.19p these formulae are symmetric with 
respect to permutations of indices and therefore all coordinates have the same 
rights for changes. Thus, the state with 

n n n 

pr = i/nj2Pk, /?r = i/^E^'^' Qr = i/^E^'^' ^ = 1--^ (3-20) 

k=l k=l k=l 
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Figure 2: The equilibrium fixed point in the phase space {Pf ,Ri) 
is fixed. □ 

Surely the above Theorem 13.11 has the computer illustration. If in the concrete 
model one put the starting coordinates Pi,Ri, Qt of vectors P, R, Q equal to 
middle-arithmetic value P-*^^, R^^, QT'i then they do not change for any > 1. 
Thus, the corresponding state of the system is equilibrium. 

3.2 The equilibrium state is stable 

The fixed point from Theorem 13. II is in fact a one-point attractor. This result we 
formulate as follows 

Theorem 3.2 The equilibrium state of the conflict triad dynamical system Ii3.16\) 
which is determined by the arithmetic mean values of the coordinates P^^, R^^, Q°^, i = 
l...,n (see Ii3.20\) ) is a one-point attractor. 

Here we do not cite the formal mathematical proof of this theorem, but only 
the computer illustration. 

Example 1. Consider the case with four regions, i.e., the number of coor- 
dinates n = 4. Let us choose the following meanings of parameters in formulas 
dSHD, (M): 

d^ = d3 = 0.09, d2 = 0.01, a = 0.1, b = 0.6, c = 0.1. 



15 



5500 6000 6500 7000 7500 8000 8500 9000 9500 



Figure 3: The cyclic attractor in the phase space {Pf, Qi) 

Note that these values are guessed "by hands". Under their unreasonable replace- 
ment, the model can "fall into pieces" (in particular, some coordinates will go to 
infinite or become negative). Further, let the population of biological species at 
initial moment of time has the following quantitative distribution along regions: 

Pi = 9000; P2 = 5000; P3 = 2000; P4 = 12000. 

We put the distributions of the vital resource and the infection density to be 
essentially smaller: 

Pi = 30; P2 = 80; P3 = 50; P4 = 10; Qi = 5; Q2 = l] Q3 = 2; Q^ = 4. 

The direct analysis of this computer model confirms Theorem 13.21 We test the 
behavior of all coordinates to = 2000 step. We observe, they demonstrate the 
following features. Each coordinate of all three vectors goes to an appropriate 
value when N — )■ 2000. In the phase spaces we get the perfect spirals which twist 
to the equilibrium point: P^^, P^^, Q^*^, i = 1..., n with the arithmetic mean values 
of the initial dates. 

It is important, that in the phase spaces {Pf^ , Rf) and {Pf^ , Qf) we get 
different spirals which twist in the opposite directions to points (7000,42.5) and 
(7000,3), respectively. We remark that values P^^, Pf^ and Qf oscillate for a 
long time. They simultaneously approach under N — >■ 2000 to the fixed means 
Pj°° = 7000, P°° = 42.5 and Qf = 3 which exactly are the middle values by 
regions of the starting distributions of the proper substances (see Fig. 1-2). 
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Figure 4: The almost square attractor in the phase space {P-i^ , Q2) 

Considering of numerous examples, in particular, varying the values of the 
initial coordinates demonstrate the stable character of the equilibrium state. No 
doubt, the state of the conflict triad defined by the middle-arithmetic values is 
attractive for all trajectories close to this state. Thus, the proper point is a local 
attractor. Moreover, even large enough change of coordinates for the vector of 
biological populations (for example, the replacement Pi = 9000 by Pi = 900) 
does not destroy the attracting property to the equilibrium state. 

3.3 Cyclic attractors 

The equilibrium state from the previous example (the locally stable fixed point) 
is not a global attractor. In particular, this state is vanished under the large 
enough change of single coordinate of the vector that corresponds to an epidemic 
infection. Of course, under any small variations of coordinates, the system does 
not leave the attracting phase to the equilibrium state. Nevertheless, the replace- 
ment Q2 = 1 by (^2 = 4 transforms our system into another behavior phase. We 
observed, in particular, the phase of attraction to the cyclic orbit. Under the 
change mentioned above all coordinates are not attracted to a fixed point (the 
equilibrium state) but approximate a cyclic trajectory. For more details, let us 
analyze the following example. 

Example 2. The loss of the equilibrium state (a stable fixed point) may take 
place under some changing of the initial coordinates and even only one of them. 
In particular, even changing of the initial values of any substances. So, replac- 
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Figure 5: The periodic oscillation of the virus coordinate Q2 

ing R2 = 80 by R2 = 40 we obtain the appearance of the cychc attractor in the 
phase spaces {Pl^,Rf) and {Pl^ ,Qf) (see Fig. 3). That is, they are achieved 
quickly enough, aheady at = 400 step. It is interesting that cychc oscillations 
of Pj^, and take place around points which are shifted with respect to the 
arithmetic mean values of the initial coordinates (see the previous example). 

How to explain this shift? There are also additional questions. For example, 
why a decrease of full mass of vital resource (from 170 to 130) leads also to 
passing of the system into the new phase? It is not attracted already by the fixed 
point. The evolution trajectory of the conflict triad approximates the cyclic orbit. 
They are attracted rather quickly to a cycle of the egg-like form. That is, in the 
limit all trajectories oscillate around fixed points which are shifted in comparison 
with the initial mean values. All these questions are open problems, although in 
[2] we made attempts to give some interpretations to the phenomena mentioned 
above. 

Example 3. The exponential increasing of the initial coordinates Qi causes 
the appearance of almost square attractor in the phase space. That is, Qi oscil- 
lates from zero to the maximal value, nearly 120. In particular, the replacement 
Q2 = 10 by Q2 = 100 induced the oscillations of large amplitude, up to 200! (See 
Fig. 4). In the same time the coordinates of resource vector Rf have periodic 
oscillations without achievement of the absolute maximum. 

These observations confirm an interesting practical effect. The substance 
which corresponds to negative threats (an epidemic infection) has large influence 
to the behavior of biological population. So, the relatively small increase of 
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Figure 6: The periodic oscillation of the resource coordinate Q2 

initial values of coordinates Qi substantially multiplies the negative effect on other 
substances. Thus, the increasing of threats presses not only to the existence of 
biological population, but also at the resource environment. In particular, it 
makes impossible the achievement of maximal values by them (See Fig. 5,6). 

3.4 The wave of cyclic attractors 

This phenomenon arises up under a certain increasing of parameters di,d2id^. 
We represent it at Fig. 7 where few cycles are observed and they imposed one 
at another. These cycles have different periods. In this case the trajectories are 
approximated consecutively to one of cyclic attractors, but only for some time. 
We call this picture the wave of cyclic attractors. Note that similar behavior of 
complex biological dynamical systems appear in practical situations, for example, 
in a case with several epidemic sources. 

The cyclic attractors for orbits can have not only egg-like shape, but also 
considerably more difficult geometrical structures (see Fig. 8). In particular, the 
cyclic attractor similar to the Carno-cycle (see Fig. 9) describes the behavior of 
a pair " virus- resource" . It has close analogy with work of a move aggregate that 
consumes a certain resource, but survives resistance and have to come back at 
the starting position. 
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3.5 The quasi-chaotic behavior 

At the computer model that corresponds to Fig. 8 we observe an obvious evolution 
non-balanced for biological population and infection, we call this dynamical phase 
as the quasi-chaotic behavior. Thus, any regularity is absent here and one can 
not find even a slightly noticeable low in the behavior. We may find some analogy 
with a chronic hidden disease, that is healed, but not cured. It increases non- 
periodically but does not reach a critical stage in an organism which passed some 
threshold in its development. In this situation it is impossible to go back to the 
state of attracting to a stable point of equilibrium without external influencing 
(See Fig. 9). 

More deep analysis shows the presence of bifurcation points in some zones of 
phase space. Depending on parameter meanings di,d2,ds and a,b,c and initial 
values of coordinates, the evolution of the whole system passing this point may 
sharply change its direction. In particular, it is not structurally stable and can 
already not reach the equilibrium state. 
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Figure 7: The multi- dimensional attractor in the phase-space {P^ ■, Q^) {di = 
0.001; ^2 = 0.000001; 4 = 0.0012; a = 0.1; 6 = 0.6; c = 0.1; Pi = 9000; Pa = 
5000; P3 = 2000; P4 = 5; Qi = 50; Q2 = 2; Q3 = 1;Q^ = 40; Pi = 30; P2 = 
40; P3 = 50; P4 = 10) 



4 Discussion 

The main result of our work consists in the construction of complex models which 
demonstrate the presence of all basic phases for coexistence of typical conflict 
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Figure 8: A Carno-type cycle in the phase-space {with = 0.0017) 

triad substances: the state of dynamical equilibrium (a stable fixed point), the 
phase of the limit cyclic attractors (trajectories with periodically pulsating evo- 
lution for each of substances), the area of bifurcation points, that corresponds, in 
particular, to series of the extremely complex behaviors with transitions between 
some cascade of cyclic attractors, and finally, the phase close to chaotic which 
we call the quasi-chaotic behavior. All these phases for coexistence of conflict 
substances one can meet rather often in reality. So, the state of dynamical equi- 
librium is typical for the simultaneous coexistence of different kinds of bacteria 
inside a living creature, when the highly dangerous microbes exist in a healthy 
organism but do not exceed a critical concentration and cause obvious rejections 
(that is an ordinary dynamical equilibrium). Although the considerable enough 
external impact is able to violate such equilibrium state and to result in creation 
of priorities to some kind of bacteria and symptoms of a certain illness. Other typ- 
ical phase, the cyclic attractors, observed, for example, when the season (cyclic) 
flu epidemics are caused by exhausting of the valuable vital resources and by a 
too much concentration (an overpopulation) of biological species in some regions. 
This phase is character of most oscillating processes in our reality. 

Certainly, the phases and states of the dynamical systems of conflict triad 
mentioned above appear only under appropriated values of parameters and start- 
ing coordinates. In the majority of situations we observe an extremely complex, 
sometimes practically chaotic behavior of our system. This all means that the 
theory developed here is in fact only an introduction to study of the conflict triad 
by the regional approach method. Of course, it can be applied as a research 
instrument for concrete tasks in conflicts with triple opponents. 



21 



4500 5000 5500 6000 6500 7000 7500 8000 8500 9000 9500 10000 



Figure 9: Several attracting sets with overlapping basins for {P-[^ , Qi) {di = 
0.95; ^2 = ds = 0.01; a = 0.1; 6 = 0.6; c = 0.1; Pi = 9000; Pa = 5000; P3 = 
2000; P4 = 12000; Qi = 5; Q2 = llQs = 2; Q4 = 4; Pi = 30; P2 = 40; P3 = 
50; P4 = 10) 

It is important that our model in concrete settings allows to determine the 
probability of infection for considered biological species by an epidemic disease 
in each separate region of common space of existence. We recall that a key point 
of our constructions is a natural division of the whole territory of existence fl 
onto a set of regions Qi, i = l...,n with studying of the local process of conflict 
interactions. Thus, the complete picture of dynamical changes takes into account 
redistribution and migration processes between regions, both for biological species 
and for infections. In essence, a model gives not only the statistical picture 
of distributions but also the quantitative characteristics for each of interacting 
substances in regions at discrete time. It is possible to provide the prediction and 
information about the power of infection risk, periods of relative safety and picks 
of sharp growth of the infection density (epidemic) using the model. Besides, it is 
also possible to pick regions with relative stability or opposite ones with biggest 
growth of disease. 

Finally we remark that it is possible to watch for sizes of relations such as 
Rf / Pf^ , Qf / Pf^ which make sense of distribution densities of positive vital re- 
sources and threats of infection for existence of biological species in region. 
These sizes are important at the decision making problem for safe of existence of 
biological species in a fixed region. The error of uncertainty of such decisions de- 
pends on the range of changes for parameters di,d2,d3,a,b,c which are controlled 
by external factors with respect to the system. 
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We hope that the dynamical model of conflict interaction between triple in- 
trinsic elements offered here may be used as a flexible tool in the problem of 
practical forecasts. Among important dynamical parameters one can take such 
ones: the speed of reproduction and spread of infection, their density and the 
local concentration, the rates of migration between regions and so on. 
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